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We propose a modification of the Hybrid- Monte-Carlo algorithm that allows for a larger step-size of the integra- 
tion scheme at constant acceptance rate. The key ingredient is the splitting of the pseudo-fermion action into two 
parts. We test our proposal at the example of the two-dimensional lattice Schwinger model and four-dimensional 
lattice QCD with two degenerate flavours of Wilson-fermions. 



1. INTRODUCTION 

In large scale simulations of lattice QCD with 
two flavours of mass-degenerate Wilson fermions, 
the mass of the fermions is still too large com- 
pared with the up and the down quark masses. 
Therefore a delicate extrapolation of the data is 
needed. Unfortunately, the numerical effort re- 
quired for the Hybrid-Monte-Carlo (HMC) [| al- 
gorithm increases as the quark mass decreases for, 
at least, three reasons the solver (BiCGstab, 
conjugate gradient) needs more iterations, the 
step-size of the integration scheme has to be re- 
duced to maintain a given acceptance rate, the 
autocorrelation times (in units of trajectories) in- 
crease. 

We propose to modify the pseudo-fermion ac- 
tion in such a way that the problem of the de- 
creasing step-size is drastically reduced. Our 
starting point is the observation that the step- 
size can be increased when the fermion matrix is 
replaced by its preconditioned counter-part Bj5|. 
This means that reducing the condition number 
of the fermion matrix allows for a larger step-size. 

Here we shall present results for the two- 
dimensional Schwinger model j(| and four- 
dimensional lattice QCD. 

2. THE SCHWINGER MODEL 

Let us start with the partition function with 
two degenerate flavours of Wilson-fermions 

Z= I D[U}exp[-S G (U)] detM(U)^M(U), (1) 
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where Sq(U) is the gauge action and M(U) = 
t — kH(U) the fermion matrix. In the HMC, the 
fermion determinant is represented by an integral 
over a so-called pseudo-fermion field (f>: 

detM+M oc J D^D^t] exp(-|AT l <p\ 2 ). (2) 

Hence the pseudo-fermion action is given by Sf = 
{M' 1 ^ 2 . The key-ingredient of our method is to 
split M into two factors, where each factor has 
a reduced condition number. For each factor we 
use a pseudo-fermion field: 

Sfi = |Af~V| 2 S F 2 = |MAf~V| 2 , (3) 

where M = 1 — kH with < k < k. 

In the HMC we have to compute the variation 
of the action with respect to the gauge-field. It 
turns out that for our modified action this can be 
done in much the same way as for the standard 
pseudo-fermion action: 

6S F 2(U, 4>) = - [y f SM X + 5M* Y] (4) 

with X = Af"V and Y = Af t_1 (<*/>+ where 
a = k/n and b = 1 — a. 

In our numerical study of the Schwinger model, 
we have applied this modification on top of even- 
odd preconditioning: M is replaced by M ee = 
i ee — K 2 H eo H oe . The modified pseudo-fermion 
action is given by 

S F 1 = \M-M 2 S F 2 = \MeeM- e 1 (}>\ 2 , (5) 

where M ee = t ee — k 2 H eo H oe . We have per- 
formed our simulations of the Schwinger model 
at the same parameters as ref. B . We have used 
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Figure 1. 

cost = (number of applications of H eo H oe per tra- 
jectory) as a function of k for L = 32, /3 = 4.0 and 
k = 0.26. The acceptance rate is fixed to m 0.8. 



the leap-frog integration-scheme. In order to sep- 
arate off effects from the gauge action, we used a 
reduced step-size (n = 4) for the gauge action |j . 
The length of the trajectory is taken randomly be- 
tween 0.5 and 1.5. We used the BiCGstab solver 
to compute M~ e lr tp and M" 1 ^. 

In all cases we have chosen (3 = 4.0. In a 
first set of runs we simulated a 32 x 32 lattice at 
k = 0.26 and various values of k. Following ref. 
Q the pseudo-meson mass is mp — 0.210(3) for 
this value of k. Each run consists of 10000 trajec- 
tories. 1000 trajectories were discarded for equi- 
libration. We have measured the value of square 
Wilson-loops up to the size 5x5 and the topolog- 
ical charge with the geometrical definition. We 
found that for fixed acceptance rate, the inte- 
grated auto-correlation times of these quantities 
do not depend on the value of k, within error- 
bars. Therefore, the performance of the algo- 
rithm can be measured just by the number of ap- 
plications of H eo H oe that is needed in average for 
one trajectory. This number is given as a func- 
tion of k in figure |l|. We see that the numerical 
cost has a shallow minimum at k w 0.22. At this 
minimum cost is reduced by a factor of about 1.66 
compared with the standard HMC (k = 0). 



Next we performed simulations of a 64 x 64 
lattice at k — 0.2570 and k = 0.2605. The 
pseudo-meson masses at these values of k are 
mp = 0.210(3) and mp = 0.124(5), respectively 
For k — 0.2605 we found a reduction of the 
cost by factor of 2 compared with the standard 
pseudo-fermion action. For k = 0.257 we see only 
a tiny improvement. This observation indicates 
that the gain in performance increases as we ap- 
proach n c . For details see ref. [||. 

3. LATTICE QCD 

We have performed some exploratory runs for 
QCD in four dimensions. Here we used a differ- 
ent set-up as in the case of the Schwinger model. 
These changes are mostly motivated by the exist- 
ing HMC-code of the ALPHA-collaboration. 

We considered the Hermitian fermion matrix 
Q = c 7 5 Af . We define Q = Q + ia 1. That 
means that the modified pseudo-fermion action 
consists of the two parts 



S F1 =^[Q 2 + a 2 t}-^, 



S 



F2 



(6) 
(7) 



In our simulations we have used 0(a) clover 
improvement and even-odd preconditioning JsJ . 

Instead of the leap-frog scheme we used a 
scheme proposed in ref. |?J that partially re- 
moves the At 2 errors. This integration scheme 
is characterised by the coefficients 1/6, 1/2, 2/3, 
1/2 and 1/6. Note that in contrast to the leap- 
frog scheme, the derivative of the fermion action 
with respect to the gauge-field has to be com- 
puted twice per elementary step. 
In a first set of runs we studied an L A system 
with Schrodinger-functional (SF) boundary con- 
ditions, where k rs k c . The parameter a of the 
modified action is set as a 1 — \J \ m in , where A m j„ 
is the smallest eigenvalue of Q 2 . Our results are 
summarised in table [l]. For L = 8 the step-size is 
enlarged by a factor of 1.4 by the modification of 
the action, while the overhead is about 1.29. For 
L = 12, we see already an increase of the step size 
by a factor of 2, while the over- head is 1.28. 

Finally we performed runs at a smaller value 
of f3 for a physically large volume. We took the 
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Table 1 

Runs with SF-boundary conditions. At is the 
step-size of the integration-scheme, N mc i the num- 
ber of steps per trajectory, "acc" the acceptance 
rate. Ncg and Ncg2 are the average number 
of iterations needed to compute (Q 2 )^ 1 ^ and 
[Q 2 + fi 2 l]"V with the conjugate gradient solver. 
The standard HMC data are taken from rcf. 
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Table 2 

Runs for (3 — 5.2, k — 0.1370 and c sw — 1.76 with 
periodic space and anti-periodic time boundary 
conditions on a 8 3 x 24 lattice. 
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N md 
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N CG 
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82% 
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10 


80% 


122 


27 



parameters for the run from rcf. ||. In rcf. || 
the leap-frog scheme was used. For At = 0.02 
and N m d = 50 elementary steps per trajectory 
they obtained an acceptance rate of 80%. 

Our results are given in table g. For these 
particular parameters, the improved integration 
scheme performs better than the leap-frog. With 
the standard pseudo-fermion action we get an ac- 
ceptance rate of 82% with At = 0.06. I.e. the 
improved scheme out-performs the leap-frog by a 
factor of 1.5. Using the modified action, the step- 
size can be further increased to 0.1, maintaining 
an acceptance rate of 80%. 

4. CONCLUSIONS 

We propose to use a modified pseudo-fermion 
action in the HMC simulation of dynamical 
Wilson- fermions. It is easy to incooperate this 
modification in an existing HMC Code. We have 
demonstrated that the modification reduces the 
numerical costs for the 2D Schwinger model as 



well as 4D QCD with clover- improvement. We 
see an improvement up to a factor of 2. However, 
we expect that for lighter quark masses this gain 
becomes even larger. 

There remain a number of open problems. We 
need longer runs for QCD to determine reliable 
autocorrelation times. We have to compare the 
modification of the pseudo-fermion action with M 
and Q in the same model. Is it useful to split the 
pseudo-fermion action in more than two parts? 
We have to study more carefully the inter-play of 
various integration schemes with the modification 
of the action. 

Finally, our idea can be combined with 
Polynomial-Hybrid-Monte-Carlo algorithm |l0| 
pT| . (See Mike Peardons plenary talk). The im- 
plementation is straight forward. The non-trivial 
question is the choice of the polynomials. 
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